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We will show how it is possible to generate entangled states out of unentangled ones on a bipartite 

system by means of dynamical boundary conditions. The auxihary system is defined by a symmetric 
but non-sclf-adjoint Haniiltonian and the space of self adjoint extensions of the bipartite system 
is studied. It is shown that only a small set of them leds to separable dynamics and they are 
characterized. Various simple examples illustrating this phenomenon are discussed, in particular 
we will analyze the hybrid system consisting on a planar quantum rotor and a spin system under a 
wide class of boundary conditions. 
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I. INTRODUCTION 

There is an increasing interest in the physics associ- 
ated to the "boundary" of a given physical system. As 
the boundary can be thought as an effective way of de- 
scribing the interaction of the system with the external 
universe, its modeling could account for a number of sig- 
nificant physical phenomena. It is impossible to cover the 
range of physics associated to boundary structures. We 
will just mention here Casimir's effect which is arguably 
one of the most conspicuous physical phenomena associ- 
ated to presence of boundaries (see for instance [As06], 
[As08] and references therein for an extensive account of 
the role of boundary conditions and vacuum structures) . 
We would like to add here the possibility of describing 
topology change as a boundary effect. This possibility 
was already considered in [Ba95] and further elaborated 
in relation with specific boundary condition in [As05], 
but it has gained new impetus because of Wilczek's et al 
[Wil2] recent contribution to it. 

In this paper we will explore how the manipulation 
of boundary conditions of composite systems allows to 
generate entangled states. More precisely, consider two 
systems A and B, and assume that the system B, which 
will be called the "bulk" or controlled system, is com- 
plete, i.e., its Hamiltonian Hb is an Hermitean (self- 
adjoint) operator on a Hilbert space Hb and its evolution 
= cxp{itHB) is unitary. However the system A, or 
"auxiliary" , is defined by a merely symmetric operator 
Ha on a Hilbert space Ha- In other words the evolu- 
tion = exp{itHA)" will not be unitary until we have 
selected (if it exists) a self-adjoint extension of the op- 
erator Ha- It is worth to point it out here that such 
situation will actually happen whenever our system A is 
defined in a bounded domain Qa in with boundary 
d^A ■ In such case the Hilbert space Ha will be the space 
£^{rt) of square integrable complex-valued functions on 
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CIa and the Hamiltonian operator will be: 



A„ 

2m ^ 



(1) 



with the Laplace-Beltrami operator defined by some 
metric r] on Qa, and V a potential function. Under such 
circumstances it c;an be shown that the self-adjoint ex- 
tensions of Ha are determined by boundary conditions 
satisfied by the functions on the corresponding domain 
[As05]. 

The main observation which is relevant for the pur- 
poses of this paper, is that if we consider the bipartite 
system defined on the Hilbert space Ha'S^Hb, the family 
of self-adjoint extensions of the symmetric Hamiltonian 
H = Ha (E) I + I <E) Hb, is much larger than the fam- 
ily of self-adjoint extensions of the standalone symmetric 
Hamiltonian Ha- As it will be discussed along the pa- 
per, many of the self-adjoint extensions of the bipartite 
system are able to generate entangled states, i.e., they 
do not define separable dynamics, meaning by that they 
will not preserve separable states. Separable dynamics 
will be characterized and it will be shown that they cor- 
respond to the trivial extension to the tensor product (in 
a sense that will be described later on) of a self-adjoint 
extension of the system A. Such sources of entangled 
states will be called boundary generated entanglement. 

It will be illustrated in simple examples how by choos- 
ing non-trivial tensor product extensions of a given 
self-adjoint extension of the system A, we obtain non- 
separable dynamics. Even more, we will show how by 
modifying the chosen self-adjoint extension, we can gen- 
erate entangled states not only between the auxiliary sys- 
tem A and system B, but even in system B itself (as long 
as it is a composite system itself). 

Such instances will be discussed first by using a toy ex- 
ample consisting of the free particle moving on the half 
line as auxilary system and a two level system as a bulk 
system. In this particular instance it will be shown that 
the ground state of the half-line (actually its only eigen- 
state) becomes entangled with the eigenstates of the bulk 
system and how such entangled state can be driven by 
modifying the boundary conditions compatible with such 
scenario. Finally, we will disciiss a "quantum compass" , 
i.e. a planar rotor possessing a spin 1/2 system sitting 
inside it. Now, two families of non-trivial boundary con- 
ditions for such system, extending in a nontrivial way 
quasiperiodic boundary conditions for the planar rotor 
[As83], will be considered and their spectral properties 
will be discussed. 



II. BOUNDARY CONDITIONS AND 
SELF-ADJOINT EXTENSIONS 

We will start by reviewing briefiy the most salient as- 
pects of the relation between self-adjoint extensions and 
boundary conditions by using the Laplace-Beltrami op- 
erator as an illuminating example. 



Given a symmetric operator T on a Hilbert space H, 
i.e., T has dense domain Vq C "H and T <Z T\ we may 
use von Neumann's theorem [Ne31] (see [Re75] for an 
exhaustive account of the theory) to describe all its self- 
adjoint extensions. Namely, we compute the deficiency 
spaces JV± = ker(T''' =f iT) = Ran(T± zl)-*- and then there 
is a one-to-one correspondence between self-adjoint ex- 
tensions of T and unitary operators K : J\f+ — )• M- . To 
any such K we associate the operator Tk with domain: 



Vk=Vq®{1 + K)N+ 



(2) 



and such that 



(3) 

In many occasions the operator T, the Hamiltonian 
of the system for instance, is defined as a diff'erential 
operator on a manifold $7 with boundary Let us 

consider in what follows, as an illustrative situation, the 
case of a free particle moving on a curved manifold O with 
Riemannian metric rj. In such cases, the Hamiltonian 
describing the geodesic motion is the Laplace-Beltrami 
operator (this is, we are in the situation of eq. (1) 
with V = Q), that in local coordinates a;*, i = 1, . . . ,n, 
n = dimf2, takes the explicit form: 



1 d 



r 



d 



(4) 



with \ri\ = dei{r]ij). It is natural to select the domain of 
this operator as C^(Int(n)), i.e., as the set of complex- 
valued functions with compact support contained in the 
interior of f2. This is a dense subspace of "H = >C^(f2), the 
space of square integrable functions with respect to the 
Riemannian volume defined by rj. A simple integration 
by parts leads to 

($,A^*) = (A^*,*) V$,4'eCr(Int(r!)), (5) 

showing that the operator is symmetric. The minimal 
closed extension of the operator A^ is defined on the do- 
main 2?o = HHP), i.e., the closure of C'^{lat{n)) with 
respect ot the Sobolev norm || • ||2, which is just the 
Sobolev space of functions of order 2 vanishing at the 
boundary together with its normal derivative. The ad- 
joint operator A!^ is the operator defined in the domain 

= {$ e C^{^) I A^$ G C'^{^)}. Such operator At 
is actually the maximal extension of A^ and certainly 
A^ C At. 

A general result on operators commuting with conjuga- 
tions shows the existence of self adoint extensions for A^, 
then the existence of unitary operators K: A/+ — >■ M- 
and the applicability of Neumann's theorem, eqs. (2), 
(3). Alternatively we may argue as follows (see for in- 
stance [As05] and references therein). Consider the re- 
striction to the boundary dO. of functions in Vq. Such 
restrictions will be denoted by </? = $ I^q. In the same 
way we define the normal derivative ip = d'^/dv |an as 
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the outbound normal derivative along the boundary. We 
will consider that both (p, (p are in L^(d^) (even though 
this is not necessarily so and care must be taken when 
analyzing particular examples) . Repeating now the com- 
putations leading to eq. (5) for elements € I?J we 
will obtain: 

($,A^*)-(A^$,'J') = (^,^)-(^,V) (6) 

where the inner product in the r.h.s. of the expres- 
sion above is the one defined in £^(90), i.e., {(p,4') = 
Jg^(fi{x)'ip{x)diigQ{x) with /ign the measure associated 
to the Riemannian metric induced at the boundary dil 
by rj. Clearly, self-adjoint extensions of A,, will be deter- 
mined by maximal subspaccs of functions $ in Vq such 
that the bilinear form given by the r.h.s. of eq. (6) will 
vanish identically for the corresponding boundary values 
(fi, (f of 

Such maximally isotropic spaces W of boundary val- 
ues can be easily characterized by computing their Cay- 
ley transform, i.e., we consider the linear isomorphism 
C: jC^idfl) e £2 (an) ^ £2 (an) e C^idn) defined by 

C{<p,^) = -^{(p-i<f,<f + i(p). 

Clearly C maps a maximally isotropic subspace W 
into the graph of a unitary operator U : £^(9fl) — ^ 
£2 (an). More explicitly {(p,(p) G W iS there exists 
U e W(£2(aO)) such that [As05]: 

if + iip = U{ip - itp) . (7) 

In this sense the space of all self-adjoint extensions of 
the Laplace Beltrami operator can be naturally identi- 
fied with the unitary group of the Hilbert space of square 
integrablc functions at the boundary of and eq. (7) 
provides the explicit description of their domains. We 
will make extensive use of this characterization in what 
follows. 



III. SELF-ADJOINT EXTENSIONS OF 
SYMMETRIC BIPARTITE SYSTEMS 

Let us consider now the case of a bipartite system AxB 
such that one of its subsystems is described by a sym- 
metric operator. In particular we will consider system A 
to be defined as in the previous section by the Laplace- 
Beltrami operator on a Riemannian manifold {flA,T]A)i 
i.e., A is a free system on a manifold with boundary. Sys- 
tem B will be defined by a self-adjoint operator Hb on 
a Hilbert space TIb with dense domain doni(iJB) = Vb- 
The Hilbert space Hab of pure states of the composite 
system is 

Hab := 'HaS'Hb = C'^{nA)®HB 

that can be identified naturally with C^{D,a','Hb): i-e. 
pure states will be considered as square integrable maps 



$ : Q.A Hb with inner product 

{^,^)ab= [ mx)Mx))ns<it^vix)- (8) 

In what follows we will use the latter identification when 
appropriate. The Hamiltonian operator will be -ff = 
— A^ 1 -|- 1 (g) Hb, that will act on states $ as: 

H^ = -Ar,^^ + HB-^, (9) 

with {Hb ■ <P)ix) = HBmx)), X e ^A. 

The natural symmetric domain T>q of the operator H 
is now Vq = Vao ^ where we are now borrowing 
the notation Vao from section II to denote the mini- 
mal closed extension of the Laplace-Beltrami operator 
defined in Qa- Again Dq can be identified in a natural 

way with Cq°{SIa] '^b)" > where the completion is taken 
with respect to the Sobolev norm of order 2. Notice that 
we can not consider the completion in the definition 
of Vq because being Vao ^nd Vb dense, it would result 
Vao®T^b = 'Ha^'Hb = 'Hab and the operator H is not 
bounded. 

The maximal extension of this operator is given by 
■^AO ® using the notation of section II again (notice 
that Hb is self-adjoint already). Computing the self- 
adjoint extensions of H is best done by using its boimdary 
data structure (i.e.. Green's formula) like in the second 
part of section 2. In fact, integrating by parts we get the 
analogue of eq. (6): 

,-Ar,,^+HB-^)-{-Ar,,^+HB-^ = - ,i>) 

(10) 

where the inner product at the boundary appearing in 
the r.h.s. of the previous equation is given simply by: 

{V,i')=j {y^{x),il}{x))nB^tJ'dr,A{x) (11) 

JdQ.A 

and if, <fi are defined as before. The space of boundary 
data is now 

c\dnA;HB) ^ cHdnA)i>HB- 

Then, repeating the argument leading to eq. (7), we 
obtain that the space of self-adjoint extensions of T, 
i.e. the space of maximally isotropic closed subspaces 
of the bilinear boundary form defined by the r.h.s. of 
eq. (10), is parametrized by unitary operators U G 
U{C'^{d^A)®'HB)- Thus, given such an unitary U, the 
domain Vu of the corresponding self-adjoint extension 
will consist of all functions $ € V\q ®iVb such that 

ip + i<p = U{tp-iip) , ip ,ip € C'^{dQA;'HB) ■ (12) 

The same result can be obtained certainly by using 
von Neuman's Theorem. Because of the decomposition 

^lo = ^Ao ® © A/'a- ) we will get: 

2?A0 ® = [T^Ao © {J^A+ © ^fA- )] ® Vb 

= {Vao ® Vb) © [{Ma+ © AT^-) Vb]{13) 
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Then the total deficiency space of the bipartite system 
will be: 

A/V ©A/l = {Ma+ ®Ma-)®'Db ^ {Na+ ®J^a-)®'Hb 
^ {Ma+ 'Hb) ® {Ma- ® Hb). (14) 

We can be more explicit by computing directly the 

deficiency spaces A/± . If wc choose an orthonormal basis 
Pn for 'Hb then any vector $ G C'^{dQ.A)®'HB has a 
unique representation as : 

^ = Y,^n®Pn. $„e/:'(ao^). (15) 

n 

Hence if p„ arc eigenvectors of Hb with corresponding 
eigenvalue A„ (we will assume for simplicity that the 
spectrum of Hb is discrete) we will get: 

(ijt =p i7)$± = (_At + HB)^<^^®PnT il^^ = 

n 

= ^(-Aj,,^^ + A„$± T i^i) ®Pn = Q 

which implies 

-At^$±T(iTAJ$±=0. (16) 
Thus must belong to the generalized deficiency spaces 
Ma,z = I H'^^+ = z^+), 

with z = (— A„ + z). However all generalized deficiency 

spaces of the form Ma,z with either Im z > or Im z < 
are isomorphic (this is, dimA/A,z is constant in both com- 
plex half-planes) and hence A/a,± = J^A,±i is isomorphic 
to A/'/i,±(iipA„) (sec for instance [Ak60]). Let us denote a 
choice for such isomorphism by : ■AfA,±(i^x„) J^A,±- 
Then we have shown that the deficiency spaces A/± of the 
operator H consist of the vectors $^ (g) /?„ with $^ € 

■^A,±(i=FA„)- The isomorphism : M± — >■ Ma,±®'Hb 
defined by 

a*(I]*^ ®Pn)-I]a^(*^)®P« (17) 

provides an explicit identification of N± with Na,±®'Hb 
as indicated in eq. (14). 



IV. SEPARABLE DYNAMICS AND 
SEPARABLE EXTENSIONS 

It is clear that if we have two complete quantum sys- 
tems A and B with Hilbert spaces of state vectors 'Ha 
and 'Hb, and Hamiltonian operators Ha and Hb re- 
spectively, then the bipartite system with Hilbert space 
vector states 'H = 'Ha®'Hb and with total Hamiltonian 
H = Ha (S) I + 1 Hb induces a unitary flow 

Ut = e**''^ = gitiHAm+IiSHB) _ git(HAi8)I)git(I®fl"B) 

= e'*"^ (g) e''-^^ = O Uf 



where , denote the individual unitary flows of the 
subsystems A and B. Then we may call a one-parameter 
family of unitary operators Uf on 'H = 'Ha'^'Hb separa- 
ble if there exist two one-parameter families of unitary 
operators and Uf on 'Ha and 'Hb respectively such 
that 

Ut = U^^ C/f . (18) 

Notice that Ut is separable if and only \iUt^ is separable 
for any separable state <I> = >I> 4 ® $5 for any t. In fact 
it is immediate to check that separable dynamics do not 
change the Schmidt index of a given state in 'Ha®'Hb- 

Now, if we are given a system H oh'Ha®'Hb which is 
obtained by means of a self-adjoint extension of the prod- 
uct of a symmetric operator on 'Ha and a self-adjoint 
operator on "Hb, can we characterize in which cases are 
we going to obtain separable dynamics? In other words, 
if [/ e U{C'^{dQ.A)®'HB) is the unitary operator defin- 
ing the self-adjoint extension eq. (12), when will it be 
separable? 

We will solve first the spectral problem for the self- 
adjoint extension of H defined by the boundary condition 
U = Ua®"^- We will assume for simplicity in what follows 
that the spectrum of Hb is discrete and nondcgencrate. 
We denote the eigenvalues and eigenvectors of Hb by A^ 

and respectively, this is Hsp^ = p^ , k = 1,2, 

An arbitrary function $ S jC^{CIa','Hb) can be written 
uniquely as: 

00 
fe=i 

hence 

00 00 

fe=l fe=l 

with<^^ = $^ \oQ^ and 93^ = d^^/du |an^. A; = 1,2, . ... 

is a unitary operator acting on 

C\dnA)®'HB, we have U^ = ^^=1 Ef=i U^'fik®Ufp^ 
and so on. In the particular instance of = f/^ 1 we 
get: 

00 

U>p = ^UA<f-^®pi 
fe=i 

and similarly for (p. Denoting as before by Hjj the self- 
adjoint extension defined by U, the spectral problem 
JJ(7$ = becomes, after some trivial computations, 
the family of spectral problems: 

H\<i>^ + Xi^t = E<ft k = l,2,... (19) 

Notice however that the boundary conditions defined by 
U become the family of boundary conditions: 

+ = UA{^t -i^t), fc = l,2,.... 
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Thus for each k we have to solve the problem: 

H\<i>^ = {E - X^)<S^^ 

Notice that if wc denote by '^f the (ugc^nfunctions of the 
self-adjoint extension of the operator Ha defined by Ua, 
this is, with the boundary conditions given in eq. (20), 
we will have: 

H{^f = Xf^f, i,f+ii,f = UA{i^f-ii>f), 1 = 1,2,... 

We will assume again, for simplicity, that the spectrum 
of the extension of Ha defined by Ua is discrete ^ . We 
denote the corresponding operator by Hu^ instead of 
Ha,Ua- We finally conclude that the spectrum of Hu 
is given by: 

E = Af + Af, k,l = 1,2,..., 

with eigenvectors ^f(^p^. Now, if $ e £^(f2A; 'Hb), we 
will have: 

k,l 

with cik = <S) Pky^)^ and if $ is separable, i> = 
^A ®^B, we will obtain: 

Cik = (*^$A)(pf,$i3> = aibk, 

with ^A = J2i (^i^f and = J2k ^kPk respectively. 
Consequently, 

e'*^^*^ $B = ^ aibke''"^ {^f ® ) 

k,l 

k,l 

which shows that the self-adjoint extension defined by 
the unitary matrix U = I/a (8) I is separable as it was 
easy to presume. 

Let us discuss now boundary conditions of the form: 

U = Ua^Ub, (21) 

with Ua e UiC^idnA)) and Ub G UIUb), i-c, separa- 
ble elements in the unitary group U{C'^{driA)®'HB)- We 
may even consider that the unitary Ub defines a sym- 
metry of the quantum system Hb, this is [Hb,Ub\ = 0. 
Now the answer is that, contrary to a simple guess, the 



^ This supposes no loss of generality for our purposes. In dimen- 
sion 1 the spectrum of any self-adjoint extension of the Laplace 
operator is discrete [We80] . In general using the spectral theorem 
[Ak60] , one can adapt this construction to the general case. 



dynamics defined by U in eq. (21) is non-separable if 

Ub^I. 

The proof of this fact is as follows. Because Ub is a 
unitary operator, it can be diagonalized and the Hilbert 
space "Hb decomposed as Hb — Wg, with Wg or- 

thogonal f/s -invariant subspaces such that: 

Ub^? ^f, y^f GWg. 

Now, because Ub commutes with Hb, Hb will leave 
the subspaces T'T's invariant too, and we will denote by 

Hb^s the restriction of Hb to Wg, s = 1,2, On the 

other hand, we have: 

oo oo 

H = jC\nA;HB) = jC\nA;^Ws) = ^jC^{nA;Ws). 

s=l s=l 

MoreovcT. the operator H]^ leaves invariant the sub- 
spaces C^{i,}A,Ws) for all s. Hence, the spectral prob- 
lem = with boundary conditions defined by eq. 
(21) is equivalent to the solution of the family of spectral 
problems: 

{H\(g)l + l^ Hb)^s = Es'S>s 

+ i'Ps = {Ua ® e'''''){ips - i'^s), s = l,2,..., 

(22) 

because $s G C?{Q.a;Ws) = C'^{VIa)®Ws; tps,(ps G 
C^{dnA)®nB, and {Ua ® Ub){<Ps - i'Ps) = {Ua ® 
e^""){iPs — i^s)- But the boundary conditions in eq. (22) 
are the same as: 

iPs+i'Ps = UA,s{V>s-i'Ps), 8 = 1,2,... 

with Ua,s = e^""UA, and in consequence the self-adjoint 
extension defined in C'^{VIa',Ws) of the Hamiltonian 
Hg — Ha + Hb,s by the boundary conditions 
Ua,s 1 is separable because of the argument in the pre- 
vious paragraph. In consequence we have obtained that: 

Finally, if $ = 4?^ (g) $b is a separable state in H, we can 
write it as: 

■oo Ns 
s=l i=l 

with G Ws- Then: 

eitHu $ = ^ ^ e''"" (23) 

s=l i=l 

^Y^Y^^itH,^^,^A^^^UH.,.^B^_ 
s=l i=l 

But now the factor e^'^' is different for each s, hence the 
eigenvalues and eigenvectors of the spectral problem in 
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C\^a;Ws), s) = 1,2,... ■ 

I ips+i^Ps=e"''UA{^s-i'Ps), s = l,2,.. 



(24) 

are different for each s. Therefore the extension Hfj^ 
is difforcnt for each s and wc; cannot factorize it out of 
the sum in the last term on the r.h.s. of eq. (23). Thus 
we conclude that if Ub I the dynamics Hu is non- 
separable. 

Actually wc can prove the following theorem. 

Theorem 1. The dynamics Hu on the product Hilbert 
space C^{VIa)®'Hb defined by the unitary operator U = 
Ua®Ub& U{/:^{dnA)§)UB) is separable iff Ub = I. 

Proof. Let assume that the dynamics defined by Hu is 
separable. Then: 

where neither Ha nor Hb have to agree with Hu^^ nor 
Hb respectively. However 

Hu = HaiSiI + IiSi Hb 

with Ha and Hb self-adjoint operators. It is also clear 
that the one-parameter group of unitary operators: 

Vt = e'*^« 

defines a symmetry group of Hu, 

[Hu,Vt] = o, yt. 

Moreover the group Vt acts unitarily in the boundary 
space £^(90a;'Hb). Then it is easy to show that neces- 
sarily: 

[U,I^Hb]=0, 

hence [Ub,Hb] = 0. But now we have a sclf-adjoint 
extension defined by a unitary matrix of the form Ua®Ub 
with \Ub,Hb\ = as in the discussion preceding this 
theorem. Then, repeating the previous arguments we will 
obtain that the dynamics is non-separable unless Ub = 
I. □ 



V. A SIMPLE EXAMPLE: THE 
HALF-LINE/HALF-SPIN BIPARTITE SYSTEM 

We will discuss now what is conceivably the simplest 
non-trivial example of a bipartite system of the kind con- 
sidered in section 3. Let the auxiliary system ^ be a 
free particle moving on the half-line K+ (O^ = R"*", 
dflA = {0}) with Hilbert space "Ha = £^(K+,(ia;) and 
the dynamics of such system is governed by the Hamilto- 
nian — ^ . The bulk system B will be a 2-level system, 
for instance a spin 1/2 system whose Hilbert space is C^. 



The dynamics is given by an arbitrary 2x2 Hermitean 
matrix Hb- We will assume that a{HB) — {Ai > A2} 
with eigenvectors a = 1,2. The corresponding bi- 
partite system A x B is defined in the Hilbert space 
H = Hi®H2 = £2(R+)gc2 ~ £2^]R+;C2). Then state 
vectors ^ gH will be written as: 



'Pi 



' P2 



$l(x)' 



(25) 

where we use the orthonormal basis {pi,/>2} to write the 
component vectors. 

As we showed before, the deficiency spaces are easy to 
compute and we get: J\f± = A/'a,± ® ~ C^, because, as 
it is easy to check, dimA/'^,± = 1, i.e., N'a,± = C How- 
ever we will work directly with boundary values which 
will prove to be more efficient. Thus, given the 
boundary values of $ will live in £^(3M"'") C^, in fact: 



and similarly 



90^ 



"$i(0)' 






$2(0)_ 







5$, 



^2 



dx 



a = 1,2 . 



x=0 



Finally: 



ip± = (p±itp 



ipi ± i^pl 
<fi2 ± 



and the self-adjoint extensions of H = — + Hb 

are characterized by unitary operators U G U{jC'^{dflA)® 
C2) ~ U{2). 

Notice that in matricial form the operator H has the 
form: 



+ Hb = 



+ Ai 







-3x^+^2 



• (26) 



We recall now that the boundary data space is given 
by £^(9]R+) (S> C^, hence according with the Thm. 1, 
separable dynamics will be given by unitary operators of 
the form U = Ua<S)1, where Ua ■ C^{dR+) £2(9R+), 
hence Ua — e*" is just multiplication by a phase. Inci- 
dentally we may recall that these are all the self-adjoint 
extensions of the system A in the half-line and they cor- 
respond to boundary conditions of the form: 



ipA + iiPA = e'^iipA - iipA) 



or equivalently 



ipA = ta,n{a/2)ipA, a ^ tt 
ipA = 0, a = tt; . 



(27) 



(28) 
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Now, because the space of self-adjoint extensions for the 
bipartite system is actually U{2) as it was shown above, 
there are many self-adjoint extensions that will define 
non-separable dynamics. 

We will consider the particular instance of self-adjoint 
extensions defined by the unitary matrices of the form: 

U = UA®y , Ua e U{C'{dW^)) , F 7^ I e C/(2) . (29) 

That despite of its form determines non-separable dy- 
namical evolution. In fact, among this class and because 
V A is just multiplication by a complex number of mod- 
ulus 1, we can just consider as the simplest non-trivial 
example a matrix V of the form 



V 



e*"i 
e'"2 



with 



(30) 



i.e. a matrix V belonging to the maximal compact tori 
inside C/(2). It is also noticeable that such V is the most 
general matrix commuting with Hb- Notice that ii ip — 
ipi® Pi+'^2® P2 e £2((9M+) (gj'Hs, then 



Uip = {l®V)ip = ipi (i 

^ Viiipi + Vi2(p2 





P2 (5 


^Vp2 


= V ■ 




= V 




}P2 





(31) 



To compute the point spectrum of the self-adjoint oper- 
ator Hjj defined by the unitary operator U = 1®V \s, 
easy. Notice that eq. (12) becomes now: 



(32) 



ipi + iipi 




e*"i 








}P2 + i'P2 




e*"2 




^2 


- iip2_ 



this is, ipa+ = e^^yJa-, a = 1, 2 or, if ai, a2 ^ tt, 

(/9a = tan(aQ/2)(^a , a ==1,2. 
Then the eigenvalue problem Hu^ — £'$ becomes 



(33) 



Hu^ = — 



' P2 



$1 ® HbPi + $2 «) HbP2 

d2 



Pl + (-d^+^2)$2' 



' P2 



E^i 'S)pi+ E^2 P2 



Then, 



Ciy/Xi — E, hence 



Moreover ipi 



dx 



\x=0 



y^Xi - E = tan(Q;i/2) or £; = Ai - tan^(ai/2) . 

Notice that E is the unique discrete eigenvalue of the op- 
erator since the rest of the spectrum {E, +00) is contin- 
uous. We can proceed similarly for the other component 
(a = 2) finding that E = \2 ~ tan2(a2/2) if -B < A2. 
In consequence, if -E < A2 we obtain the compatibility 
condition ^ 



tan2(ai/2) - tan2(a2/2) = Ai - A2 > 



(36) 



that must be satisfied for the existence of an eigenvec- 
tor with eigenvalue Ai > A2 > E. Figure 1 shows the 
space of self-adjoint extensions (ai,Q;2) with nondegen- 
erate ground state E for various values of the spectral 
gap cr = Ai — A2 of the bulk system. 

If A2 < i? < Ai, i!^ is an eigenvalue again, but this 
time the eigenvector is going to have only the a = 1 
component. We want to stress that the compatibility 
condition eq. (36) is only necessary for the existence of 
discrete spectrum. If it is not satisfied, then the problem 
has no discrete spectrum. However it still provides a well 
defined self-adjoint problem. 



0:2 



tan^(ai/2) 




FIG. 1. Phase space of self-adjoint extensions of the half- 
line/half-spin system as function of the spectral gap a pos- 
sessing a unique ground state. 



ipa = tan(aa/2)(^a 



We may start solving: 

d2$i 



a = 1,2 



(34) 



= {E-\i)^i: tpi^tan{ai/2)tpi . (35) 



Then we immediately see that if Ai < E, the solutions 
to this problem are not in £^(M+), thus Ai > i? and 



The curves defined by eq. (36) determined by the val- 
ues of (T provide families of non-separable self-adjoint 
extensions of H compatible with the structure of Hb- 
Suppose now that we select as initial state the eigen- 
state corresponding to the extension defined by ai = 



^ Remember the assumption Ai > A2 on the spectrum of the op- 
erator Hg. 
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arctany^, a2 = 0, this is $o = e~^/^^ (8) pi and con- 
sider the one-parameter family of self-adjoint extensions 
defined by the unitary matrix: 



with s, s' such that tan^ s — tan^ s' — a, this is s' = 
arctan V tan^ s — a. Then if we proceed to modifiy the 
self-adjoint extension adiabatically the eigenstate $o 
will change with s but it will remain close to the (unique) 
eigenstate of the self-adjoint extension Hjj^ , its ground 
state, that will be given by: 

= Cie-('^"''-) ^®pi-KC2e-('''"'*') ^(8)P2 , < s < 7r/2 . 

(38) 

such state $s is generically an entangled state in 
"Ha^Hb- Notice that the phase diagram of the self- 
adjoint extensions constructed in this way is periodic and 
looks as shown in Figure 2 with black dots correspond- 
ing to separable states of the form either e~^^ (8) pi or 
e"^^ 8) p2. 
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FIG. 2. Curves of self- adjoint extensions in the Abelian torus 
C (7 (2) with a single point spectrum. 



The half-line/multipartite spin 1/2 system 

We can elaborate the previous example again by con- 
sidering a system B that is already a composite system, 
i.e. 'Kb = T-Lbi'S>'Hb2 with dimHs^ = n^, a = 1,2. 
The self-adjoint operators Hb^, a = 1,2 have eigenval- 
ues A^"\ fca = 1, . . . ,nc( and a basis of eigenvectors of 
the operator Hb^ I -I- 1 8) iJsa given by 

=pi'>pi'' , (39) 



If necessary we may reparametrize s = s{t) where t is physical 
time, in such a way that < ds/dt << 1. 



where p " are eigenvectors with eigenvalues A^" . The 
eigenvalue corresponding to the eigenvector PkiM is just 
A^^^' + A^^^ Now we compute the system Ax to get 

n^nA<^nB ^c^{R+;nB,^nB2) (40) 

and we expand $ e "H as : 

J2 *fei,fc2 «)pfci,fc2 ■ (41) 

l<kc<na 

In the same way (p = J2i<k^<n^ 'PkiM ® PkiM and 
'P = Ei<fc„<„„ 'fkiM «> PkiM with (/3fci,fe2 = $(0)fe,,fe2 

and <Pfci,fe2 = ~~~5t~L-o' Filially, we notice that the 
space of self-adjoint extensions of the composite sym- 
metric operator H is given hyl4{£'^ (9IR+ ®'Hbx®'Hb2))i 
i.e. 

Mab = U{ni-n2) = U{N) , N ^ ni-n2, = dimHs„ 

(42) 

Separable dynamics correspond to J7 = C/^ x I, C/^ = 
e"'. Hence, let us choose boundary conditions leading 
to unseparable dynamics in the composite system Ax B 
and in B itself. 

We consider for instance U = lxV with V G U{'Hbi ® 
T-Lbo)- Again, choosing a simplifying hypothesis, and con- 
sidering the ordered spectrum of the Hamiltonian Hb, i.e. 
Ai = max{Ai'^) +Ai'^)} = aI^^ +Ai'^ > A2 = A^^^ + A^'^ > 

■ ■ ■ > — nrinjA^^^'' + a[.^'} we denote the correspond- 
ing eigenvectors as IIi, . . . ,11^?, then Hb^i ~ AjII;. We 
choose now the matrix V to be diagonal in this basis, 

V = diag(e*"Se'"%...,e*"") 

and repeating the computations performed in the previ- 
ous example we will get that the point spectrum of the 
operator H is given by £^ = A; — tan^ a/, / = 1, . . . , 
which imposes A'^ — 1 conditions on the parameters a/ of 
the form 

A;-tan^ Oil = A;+i-tan^ a;+i , Z = 1, . . . , iV-1 . (43) 

The previous equations (43) define a curve in the A^- 
dimensional maximal compact abelian subgroup of U{N) 
similar to those exhibited in Fig. 2. Again, a similar 
analysis as in the example of a single spin 1/2 system 
allows to conclude that an adiabatic deformation of the 
system along this curve will take a separable state, for 
instance $11 pu = $11 p^"'^'' ® p'^i^ into (maximally) 
non-separable states. 

VI. THE QUANTUM PLANAR ROTOR-SPIN 
SYSTEM 

We will consider as a final example the interesting case 
of an hybrid system that captures some properties of 
electron-nucleus systems described recently (see [Sal2]). 
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System A will be now a particle moving in the interval 
= [0,1] with measure dx, i.e. Ha = J0^{[0,l],dx). 
Unlike in the previous case, now the boundary of system 
A has two points and therefore the self-adjoint extensions 
of system A alone is going to be parametrized by matrices 
in U (2). Actually, we are going to consider a planar rotor 
with quasi-periodic boundary conditions [As83], i.e., the 
previous system with self-adjoint extensions determined 
by the unitary matrix 



Ua,, 



e'^ 



^U{C\d[Q,l])) , 



that corresponds to boundary conditions $(0) 



$(1) 



and $'(0) = e $'(1). Now we will consider as bulk sys- 
tem B a two-level system, for instance a spin 1/2 system, 
with dynamics given by Hg — fia^, where is the di- 
agonal Pauli matrix 



1 

-1 



and n is a. constant that accounts for both, the cou- 
pling constant of the magnetic field with the spin 1/2 
system and the strength of the magnetic field. Then, 
n = C'^{S\^)^C^ = /:2(5'i;C2) is the state space of 
the total system and 



H 



dx^ 



Hb 



(45) 



the total Hamiltonian. For this particular example we 
are going to change slightly the notation and will write 
the eigenstates corresponding to the Hemitean matrix 
Hb with the also usual notation for spin systems |ti)- 
Therefore i?_B I ti) — i/^lti) a-nd a particular element $ 
of the composite system T-L — 'Ha^'Hb will admit the 
decomposition $ = $''' (g) |t) + (8) As boundary 
conditions we choose U e U{<C\) (g)U{C%) of the form 
U = Ua,s ® Ub, with Ua,s as in eq. (44). 

Physically, this system can be interpreted as follows 
(see Fig. 3). One has an charged particle rotating in a 
fixed orbit [As83]. The angular momentum of this par- 
ticle depends on S. In the center of this orbit, there is a 
fixed spin that interacts with a magnetic field of strength 
fi perpendicular to the plane of the orbit. The compo- 
nent C/b of the boundary condition shall be interpreted 



as a macroscopic interaction triggered when the orbiting 
charged particle traverses an ideally infinitesimal region 
of the orbit. 

We are going now to consider two different meaningful 
situations ^ for the boundary conditions corresponding 
to subsystem B. The first situation will correspond to 
select the unitary matrix Ub diagonal in the basis of Hb , 
namely: 



(44) Ub = 



\e'" 1 

I r\ — -jVv I 



(46) 




FIG. 3. Quantum Compass 



Under this unitary matrices the boundary conditions 
defining the system take the explicit form: 

$^(0) + i$f'(0) ^ e*("+*) ($^(1) + ^$1■'(l)) 
$T(i) _ i$t'(i) ^ e^("-<5)($t(o) _ i$t'(o)) 
$^(0) + i$^'(0) = e*(-"+'5) ($^(1) + i$^'(l)) 
$^(1) - i$^'(l) = e*^-"-'^) ($^(0) - i$^'(0)) . 

One can proceed like in the previous examples and im- 
pose the above boundary conditions to the general solu- 
tion of the spectral problem eq. (45) given by 

$^(a;) = Ae''^^^"^ + Be''^'^'' 

to find the corresponding spectral function associated to 
the problem, whose zeros are the eigenvalues. By doing 
so, one obtains the following spectral function: 



* Compare them with eq. (30). 
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aa{E) = hi sillily E - fi) + 2iE sin( y^E - ^)2i^ siniy^E - fi) - - cos(J)e* 



2^sm{^/E~+Jl) ■ ■ ■ 



■ ■ + %^E - ^icos{^E - ^i) cos(Q;)e*" - 2i{E - + sin(£; - /i)e*- 

• • + 2iE sin{y/E + /i) + 2i/isin(^£; + /i) - 8^£; + /i cos((5)e"''" + %y/~E^^cos{y/~E^^) cos(a)e" 

• • - 2i(i; + + 1) sin(£; + /i)e"""^ 



Finding the zeros of this trascendental function is a task 
that has to be done numerically. However, this task can 
be challenging, especially because CTq. is very close to van- 
ish in some regions. Moreover, the information about the 
separability of the dynamical evolution depends on the 
eigenfunctions of the problem, that according to section 
IV need to admit a factorization ® pt, where the in- 
dices I and b are independent of each other. It is worth 
to remark that it is not sufficient for guarantying separa- 
bility that the eigenfunctions are separable. For all these 
reasons, in order to check that the above problem is not 
leading to separable dynamics, we will take the approach 
introduced in [Ibll]. There, an algorithm based on the 
Finite Element Method is introduced that is able to solve 
the spectral problem for any self-adjoint extension of a 
ID Schrodinger problem. Then it is enough to use the 
isomorphism £^([0, 1]) ^ ~ £^([0,1]) © £^([0,1]) to 
rewrite the problem given by eq. (45) into a form that 
can be handled by this numerical procedure. Figure 4 
shows the eigenfunctions corresponding to the 6 smallest 
energies returned by the algorithm for /j, = 10, 6 — 7r/2, 
a = tt/2. The particular values of the energies are not 
shown because they are not relevant for the discussion. 
As it can be appreciated, the eigenfunctions are separa- 




ble states in this case. However, comparing the functions 
El with £"3, E2 with E4 and E^ with Eq^ it can be seen 
that the imaginary parts (plotted in red) are changing 
sign between them. Hence we conclude that we have 
non-separable dynamics. 

Now we consider a different situation where the unitary 
matrix Ub is taken antidiagonal with respect to the given 
basis of Hb and given by 



Ub = 







-^P 



(47) 



In this case the boundary conditions defining the system 
take the form: 



$f(0) + i$^'(0) = e'('3+'5) ($^(1) + i$^'(l)) 
$^(1) - z$^'(l) = e'('3-*)($^(0) - i$^'(0)) 
$^(0) -h i$^'(0) = e'(-'5+*) ($^(1) + i$^'(l)) 
_ ^ e''^-f}-s) ($t(o) _ i$t'(o)) . 

Again, one can compute the spectral function associated 
to this problem and we get: 



0.25 0.5 0,75 1 0.25 0.5 0.75 1 0.25 0.5 0.75 1 



(Jf}{E) cx V-B^ - cos{^/E - cos{\/E + p) 

E sin{^E- ii) sm{^E + ^i) ■ ■ ■ (48) 

yjE - fi^jE + fj.cos{25) . 

Surprisingly, the spectral function does not depend on 
the parameter /3 in this case, but the eigenfunctions do 
so. In Fig. 5 are plotted the eigenfunctions corresponding 
to the case /i = 10, 5 = 7r/2, /? = tt/2. One can easily 
appreciate that they are non-separable. 





0.25 0.5 0.75 1 0.25 0.5 0.75 1 0.25 0.5 0.75 1 



FIG. 4. Eigenfunctions of the 6 lowest energy levels for the 
case fi — W, S — tt/2, a = tt/2. In each graph "1>^(3;) and 
$^(a;) are plotted simultaneously. Real parts are plotted in 
blue, imaginary parts in red. 



VII. CONCLUSIONS AND DISCUSSION 

Along the paper we have shown that it is possible by 
manipulating boundary conditions of a class of bipartite 
systems, transform a separable state into an entangled 
one. The reason for this phenomenon lies in the exis- 
tence of many self-adjoint extensions of a bipartite sym- 
metric system that lead to non-separable dynamics. We 
have been able to characterize all boundary conditions 
leading to separable dynamics in a class of symmetric bi- 
partite systems and we have discussed a number of simple 
examples were this phenomenon is explicitly shown. 
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El E2 E, 




0.25 0.5 0.75 1 0.25 0.5 0.75 1 0.25 0.5 0.75 1 



E4 Es Es 




0.25 0.5 0.75 1 0.25 0.5 0.75 1 0.25 0.5 0.75 1 



FIG. 5. Eigenfunctions of the 6 lowest energy levels for the 
case fi = 10, 5 = iv/2, /3 = 7r/2. In each graph "I>^(a::) and 
are plotted simultaneously. Real parts are plotted in 
blue. The imaginary parts vanish identically in this case. 



The systems exhibited are hybrid systems and one of 
them, the control or auxilary system, must be symmetric 
but not self-adjoint or Hermitean. The most remarkable 
fact about this class of systems is that the space of self- 
adjoint extensions is much larger than the space of ex- 
tensions of the control system and it incorporates bound- 
ary data that affect simultaneously the control and the 
controlled or bulk system. The controlled system above 
has a unitary dynamics, but together with the control 
system it becomes non-separable, hence taking the par- 
tial trace wit respect to system A, will not give us back 
the original dynamics . These ideas can be used to 
generate entangled states in a precise way, or to help to 
preserve entanglement without actually interacting with 
the "bulk" of the controlled system. The relation of these 
ideas with recent work on adiabatic computation and ro- 
bust entanglement in hybrid systems will be pursued in 
the future. 
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